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I. INTRODUCTION 



The Skyrme model is considered as the low energy limit of the quantum chromodynam- 
ics(QCD), it models QCD in the classical or large number of colors (Nc) limit and baryon 
is regarded as the soliton in the pion field HQS , Upon quantizing a slowly rotating 



Skyrmion, static property of nucleons and A have been calculated with the results in agree- 



ment with experimental data within 30% 



m 



Recently it has been widely used to discuss 



the exotic hadrons [s], ^, 0|- The minimal version of the model consists of the following 
Lagrangian terms: the non-linear Sigma term with chiral order O(p^) and the Skyrme term 
with O(p^). Even though the minimal version of Skyrme model (Min-SKM) can be regarded 
as a successful phenomenological model in spite of its simplicity, it can not be used to study 



the problem of quark spin contents of proton or EMC effects [8|, |9|, Il0| which are important 
QCD effects in baryon physics. This is a very unsatisfied defect for Min-SKM. In order 
to cure this disease, additional terms with 0{p^) or high orders have to be added into the 
model's Lagrangian to construct modified Skyrme Models. Among them, the simplest one 



is the model with the Min-SKM Lagrangian plus only one additional BnB^^ term 



9|, where 



Bfj^ is the baryon current ( or Goldstone-Wilczek current ). Hereafter, we shortly call this 
simplest modified Skyrme model as Mod-SKM. It is expected that the Mod-SKM should be 
more realistic than Min-SKM. To discuss this issue, and to fix the parameters in Mod-SKM 
is one of the aims of this paper. Moreover the Mod-SKM can be obtained by considering 
;he infinite uo mass limit of the vector meson uj term of the chiral Lagrangian studied in Ref 



An interesting application of the Skyrme model is the investigation of t 

3,Q, 



le baryon-baryon 



interaction, especially the nucleon-nucleon(A^A^) interaction []J ill3|, ^,bL6|,bJ. The 
Skyrme picture gives us a qualitative understanding of the principal features of the A^A^ 
interaction: it has the correct long-range one pion exchange potential which dominates the 
tensor force. There is a strong short range repulsion, and finally there is a pronounced 
central attraction at intermediate range, albeit weakly attractive while comparing with the 
phenomenological potential. However, the recent development of obtaining the A^A^ interac- 
tion from the Skyrme model has shown that the combined effect of the careful treatment of 
the nonlinear equations and the configuration mixing is to give substantial central midrange 
attraction for the A^A^ system that is in qualitative agreement with the data 16|. 
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The NN and A^A^ potentials have been investigated by means of the Min-SKM and 
the algebraic methods in Ref. fl8 . 



19 



20|. The phenomenon and puzzles in the baryon- 



antibaryon physics have attracted much attention recently due to the remark 



baryon-antibaryon enhancements in the J/ ip and B decays 21 



22 



kab 
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25 



e discovery of 



26|. The NN 



interaction and the possible nua 
the constituent quark model 127 



eon-antinucleon bound states have been investigated from 



28 



29l |. In the Skyrme model, the interactions between 
classical Skyrmion and antiskyrmion, i.e., SS, were explored in jy, Q]- In the present paper, 
we shall study the NN potential using the Mod-SKM and following the methods developed 
inRefs. 16, 18,Q. 

It is well-known that phenomenologically the A^A^ potential is not as well established 
as the A^A^ potential. At distance less than about 1 fermi the interaction is dominated by 
annihilation. However, at larger distance, a meaningful potential can be defined and studied 
either by G parity transformation on the A^A^ meson exchange potential or phenomenologi- 
cally. 

We will compare our Mod-SKM's results to some phenomenological potentials. The 



, l3Q 



3l|. We will 



B^B^ term in Mod-SKM reflects the effect of u meson exchange pjj, 
see that at large distance, where the product ansatz makes the best sense, the potentials 
based on the Skyrme model agree qualitatively and, in most cases, quantitatively with the 
phenomenological interactions. At intermediate and short distance, we do less well, but at 
these distances the product ansatz is not valid. However, the results are still suggestive. 

In the following section, we give the Mod-SKM's Lagrangian, then reproduce a number 
of static properties of single baryon which is both qualitatively appealing and quantitatively 
satisfactory. In Sec. HI we study the skyrmion-antiskyrmion interaction in Mod-SKM, and 
project them to the nucleon space by the algebraic methods [l8, [l9, [2^. We also consider 
the effects of rotational excitations by including intermediate states A and A, and evaluate 
the corrections to the A^A^ potential in perturbation theory. Sec. IV closes this paper with 
some discussions related to the present study. 



3 



II. THE MODIFIED SKYRME MODEL AND THE STATIC PROPERTIES OF 
SINGLE BARYON 



The Skyrme model lagrangian is generalized to include additional B^B^ term which sim- 
ulates the effects of uj meson, and this modified Skyrme model lagrangian provides a better 
description of both the single baryon static properties and the low energy A^A^ interaction 



301 . l3ll ]. This lagrangian has the following form 



£ = ^T,{d,Ud^U^) + J-^Tr([t/t9,f/, U^dM') + \mlFlTT{U - 1) - ^^^^^^^^m (1) 

where [/ is a SU{2) valued field U = exp[2iTa'n'a/ F'n], B^ is the topological current B^ = 
j^e^'^°''^TT[{Wd^U){WdaU){Wd/3U)], and e, F^, m are parameters to be determined. The 
first term is the lagrangian of meson fields in the nonlinear sigma model and the second term 
is the so called Skyrme term which stabilizes the soliton. The third term is the pion mass 
term and the fourth term is the additional B^B^ term. U transforms as U ^ U' = LUR^ 

T, X SU{2)fi, where both L and R are SU{2) matrices. 



under the chiral group SU{2\ 

The chiral soliton model ll| where, as an alternative to the Skyrme term, the vector 
meson uj term (3uj^B^ stabilizes the soliton, provides a support for the interpretation of the 
Bfj_B^ term which emerges in the limit ^ oo. In traditional nuclear interaction theories 
within the potential framework which are based on the single meson exchange mainly, it is 
shown that A^A^ system is more attractive than the A^A^ system due to the fact that in the 
theories there is a strong uj exchange, so inclusion of this term which models the effect of 
the UJ meson could help in a better description of the NN system. Furthermore, the study 
of the quark spin content also support that we should add this six derivative term in order 



to yield a spin content consistent with the present experiment (8|, |32]. Generally terms in C 
with more than two time derivatives can lead to pathological runaway solutions when the 
adiabatic approximation is relaxed and present obvious difficulties in quantizing the theory. 
But the Lagrangian of Mod-SKM have, at most, two time derivatives, hence there is no such 
difficulty. 

For the case with single static Skyrmion, we use the so called hedgehog ansatz: 

f/o(r) = exp[2r,f,F(r)] (2) 

where F{r) is the chiral angle which minimizes the static soliton energy subject to the 
boundary condition -F(O) = tt and F{oo) = 0. From Eq. ([TJand Eq. ([2]), we get the mass 



of classical soliton: 



lE^ r{x^F'^ + 2S^ + 4S\2F'^ + ^)+ 2/iV(l - C) + u^^F"} (3) 
2 e Jo 



with 



z/' = ^^, 5 = sinF, C = cosF. (4) 

In Eq. ([3]), the term proportional to u"^ comes from the B^^B^ term and is absent in the 
conventional Skyrme model. Minimizing Ms with respect to F, 5Ms = 0, we have the 
following equation for F, 

From the above equation, we can see that the chiral angle F asymptotically tends to the 
following expression when r goes to infinity, 

The coefficient A is related to the pion-nucleon coupling constant Qt^nn through 

g-KNN — — \ I ) 

Associated with the chiral symmetry, the vector current Jy^ and axial vector current J^" 
can be obtained from the Skyrme lagrangian Eq. ([1]) following the standard procedure, 

8 2 8e 2 

Q AT" 

+ [^,dM^UWU^U,d''U^U]} - ^ e'^'^^^B^Til^idaUU^df^UW - d^WUd^U^Um 
8 2 8e 2 



Q 1\T 

-[^,dM^U][d^WU,d''U^U]} - e^''''^B,TT[^{d^UU^dpUU^ + d^U^UdpU^U)l9) 



The classical field configuration of the hedgehog form does not have definite spin and isospin. 
However, nucleons carry both spin and isospin, and in any reasonable model of nucleons the 
appropriate spin and isospin states must appear. Following the conventional way, we perform 
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the collective coordinate quantization. We make a time dependent SU{2) rotation of our 
static soliton solution, 

U{x) = A{t)Uo{x)A\t) (10) 

then 

L = -Ms + / Tr(9oAt(t)9oA(t)) (11) 

where A{t) G S'f/(2)-matrix is the collective coordinate, and / is the moment of inertia, 
which is given by 

/ = r dx{S^[x^ + A{x^F'^ + 5^)] + v^S^F'''} (12) 

F^e-^ 3 Jo 

If the SU{2) matrix A{t) is parameterized by A{t) = + iTna^ with Og + Ylli=i = ^i ^^e 
Hamiltonian is 

^^-M.-ii:(A)^^M,4^M,4 (13) 

Noting that the I in the denominator is the moment of inertia, S and I are the spin and 
isospin operators respectively. As in Ref . 3| , we can calculate the static properties of single 
baryon. In going from the classical results to the quantum results for rotation operators we 
must symmetrize them 321], i.e., we perform the Weyl order of these operators. 
From Eq. f|T3l) the masses of nucleon and A respectively are 

M;v = M, + |- , Ma = M, + ^ (14) 

The isoscalar and isovector mean square electric radii are 

The corresponding proton and neutron mean square charge radii are 

{r')E,, = \{{r')E,=o + {r')E,i=i). 

{r')E,n = \{{r')E,i=,'{r^)E,i=i) (17) 

After some somewhat tedious but straightforward calculations, we can obtain the proton 
and neutron's magnetic moment which are respectively, 

/^p = 2M^(^(r2)^,,=o + ^) , 

/^n = 2M^.(^(rV/=o-^) (18) 



In the above, we have symmetrized the rotation operators, and the proton and neutron's 
magnetic momentum are defined through 

1 



(p,l/2|/i3|p,l/2): 

(n, l/2|/i3|n, 1/2) 



2Mn 
1 



2M, 



N 



-/in 



(19) 



After some lengthy calculations, we can also get the axial coupling constant 



9A 



1^ rrfa;x^{^[l + 4(F- + ^)]+F'(l + ^)} 
9e^ JO X x^ 

g2 



-2CS 



8S^ 



12 



POO OC^ 

/ dxS^F'{{eF^fr-^F' 
Jo 



3x 



3x2- 



{AxCSF' + 



18/2 



} 



(20) 



There are three parameters in the modified Skyrme model, i.e., e, Fj^, m, the pion mass 
is rrij^ = 138MeV, and the number of color Nc = 3. In the conventional Skyrme model, 
there is always a confiict between e- and F^-datum input setting for giving correct nucleon 

But a satisfactorily 



and A masses or giving correct strength of the pion tail [l4, Q, fisl - 
simultaneous description of the nucleon, A mass and the strength of the pion tail is possible 
by properly choosing the parameters CjF^^, m in Mod-SKM. Throughout our calculation 
we choose the three parameters as e = 19.48, = 129.11MeV, m = 420MeV, this 
parameter setting gives Qt^nn ~ 13.5 through Eq. ([7]), which leads to the correct one-pion 
exchange potential of A^A^ interaction as the distance tends to infinity. The connection 
between the Mod-SKM and the chiral soliton model including uj meson 11[, allows us relate 
the parameter m to the coupling /5, i.e., m = J | , and the best fit of the parameters 



in Ref. 11| gives m ^ 298.8MeV, which is not too far from the value of m in this work. 
The static properties of single baryon are summarized in the Table I, and the results of 
conventional Skyrme model ^ as well as the experimental values are also given in Table I. 

In Table I we can see that the prediction of the Mod-SKM are closer to the experimen- 
tal values than those of the conventional Skyrme model ^, so we expect that Mod-SKM 
provides a better description to other static properties of baryons including the low energy 
A^A^ interaction. 
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TABLE I: Static properties of single baryon in the modified Skyrme model compared to those in 
the conventional one and to the experimental results. 



Physical Quantity 


Modified Skyrme Model 


Conventional Skyrme Model 


Experiment Results 


Mn 


938.9 MeV(input) 


938.9 MeV(input) 


938.9 MeV 


Ma 


1232 MeV(input) 


1232 rmMeV(input) 


1232 MeV 


AU 


138 MeV(input) 


138 MeV(input) 


138 MeV 


e 


19.48 


4.84 






129.11 MeV 


108 MeV 


186 MeV 


/^2a/2 


0.71 fm 


0.68 fm 


0.72 fm 


/^2a/2 

\^ lE,I=l 


1.04 fm 


1.04 fm 


0.88 fm 




2.01 


1.97 


2.79 




-1.20 


-1.23 


-1.91 


9A 


0.82 


0.65 


1.24 



III. ADIABATIC NN INTERACTION 
A. FORMULATION 

We now study, in the product ansatz, the interaction energy of the Skyrmion- 
antiSkyrmion(S'S') system, which is a function of separation between S and S and the 
relative orientation. This interaction energy can be calculated numerically. We rotate the 
two solitons independently in SU (2) space, 

Uo{r - R/2) ^ AUo{r - R/2)A^ , 

[/J(r + R/2) ^ BU^{r + K/2)B^ (21) 

where both A and B are SU{2) matrices. In order to obtain the static A^A^ interaction, 
we describe the A^A^ configuration with the product ansatz (exact in the large R limit) as 
follows, 

U{r) = AUoir - K/2)A^BU^{r + R/2)B\ (22) 

where one baryon located at R/2 and the antibaryon at —R/2. Retaining only the potential 
energy density in the modified Skyrme lagrangian ([1]), the energy in the field of Eq. fl22|) is 

8 



the same as in 

U{r) = Uo{r - K/2)CU^o{r + R/2)C^ (23) 

where C = A^B = C4 + ir ■ c is a SU (2) matrix too. Discarding non-static terms containing 
time derivatives, the static A^A^ potential is defined by, 

V(R, C) = - [ d^xC[U{r)] - 2Ms (24) 



V(R, C) can be written in the notation of Vinh Mau et al. 13|] as 
V{Yi,C) = Fi(i?) + \/2(i?)c^ + ^3(i?)(c-R)' + \4(i?)4 + V^5(^)c^(c-R)' + F6(^)(c-R)' (25) 

where Vi{i = 1 — 6) are functions of R. Generally, for SS, the symmetry R — > — R is broken 
by the product ansatz, and we need three additional terms for a consistent expansion, 

F(R, C) = ViiR) + V2iR)cl + VsiR)ic ■ + V^iR)^ + V5{R)cl{c ■ R)^ + Ve{R){c ■ R)^ 
+V7{R)c4{c ■ R) + Vs{R)cl{c ■ R) + Vg{R)ci{c ■ Kf (26) 

These terms odd in R are artifacts of the symmetry of the product ansatz and should be 
discarded. One can use the symmetrized energy ^(^'^)+^(-'^''^) extract Vi(-R) to Vq{R), 
since the VjIR) to Vg{R) terms drop out in this combination. 

Next, we have to map the Skyrmion-antiSkyrmion(S'S') interaction to the nucleon- 
antinicleon(A^A^) interaction. This problem has been tackled in various ways by various 



groups for the A^A^ case 



13 



18|. Each o 



be cast in the form of the algebraic model 



the forms used in these works can always 
. So we will also use the algebraic method 



for mapping the SS interaction to the A^A^ interaction 



191]. This method allows 



us to study both the large A'c limit, as well as to include the finite A'c effects expli citly 



19|, 



in a systematic way. Most of the formulas given below can be found in Refs. 16 
however for the sake of completeness we remind here the important ones. 

The algebraic model consists of two sets of U{4) algebras, one for each Skyrmion (or 
antiSkyrmion), as well as the radial coordinate R. This method was developed in Ref. 



16 



ISj for the A^A^ system, and also generalized to the A^A^ system 19 



20|. In large A^c 



limit, the SS interaction can be expanded in terms of three operators: the identity, the 
operator W and the operator Z, 



Z = T^;r^^[?,R,R,-5,,]/Nl 
9 



(27) 



Here a and /3 label two different sets of bosons, used to realize the f/(4) algebra, and T is 
an one-body operator with spin and isospin 1. The semiclassical (large Nc) limit of these 



operators can be given in terms of R and C = C4 + ir ■ c as jl8| 

WM.B) = 3cl-c^ = 4cl-l, 
Z^i{A, B, R) = 6(c ■ R)2 - 2c^ = 2c^ - 2 + 6(c ■ R)^ (28) 

The SS interaction can be expressed as 

V{R, C) = vi{R) + V2{R)W,i + v;{R)Z,i + v^{R)Wl + v,,{R)WaZ,i + v^{R)Zl (29) 

in the semiclassical limit. We can obtain the relations between Vi and Vi{i = 1 — 6) by 
comparing Eq. fl25p and Eq. 



Vi{R) = vi{R) - V2{R) - 2vs{R) + v^{R) + 2^5(7?) + 4ve{R) , 
V2{R) = 4v2{R) + 2vs{R) - 8v^{R) - lOv^iR) - 8ve{R) , 
VsiR) = QvsiR) - 6v,iR) - 2Ave{R) , 
V^{R) = IGv^iR) + 8v5{R) + 4vq{R) , 
V^iR) = 2Av5iR) + 2Ave{R) , 

V<,{R) = 3Qv^{R) (30) 

Six independent choices of the matrix C can yield enough independent linear equations to 
determine Vi{R){i = 1 — 6) or equivalently Vi{R){i = 1 — 6) through Eq. (!30|) . and the 
numerical results for Vi{R){i = 1 — 6) coming from the Mod-SKM are shown in Fig. 1. From 
this figure it can be seen that the first three term vi{R), V2{R) and v^lR) are dominant. 
It seems a good approximation to neglect the interaction terms which are nonlinear in the 
expansion of operator W and Z. In the following discussion, we will mostly concentrate on 
the first three terms, then the leading term in this expansion is given by the following form: 

V(R,C) =vi{R)+V2{R)W + V3{R)Z (31) 

The algebraic operators W and Z have simple expectation values for the nucleons 3] 

Nr 

{N\T^,\N) = -^P^{N\T^am) , 

{NN\W\NN) = -P^{NN\(T^ ■ a^T^ ■ t^\NN) , 
9 

{NN\Z\NN) = Ip^{NN\{3(t'^ -ilcr^ -R-cr^ ■cr^)T^ ■t'^\NN) (32) 
9 
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Here P/v is the finite Nc correction factor P/v = ^ + -By using Eq. (132]) we take tlie NN 
matrix element of the interaction and evaluate the A^A^ potential, which only contains three 
independent multipole component, i.e., the central part Vc, the spin-spin part Vs, and the 
tensor term Vt. 

V^(o)(R) = Vc{R) + Vs{R)(T^ ■ a-'^T^ ■ + Vt{R){?>CT^ ■ Rcr^ . R - cr^ ■ cr2)r^ ■ (33) 
with 

V. = v,. K = ^, ^. = ^ (34) 

The NN potential in the above is calculated by projecting Eq. fl3T]) to the nucleon degrees 
of freedom only, and this is the correct procedure for large separation. However, at short 
distance the nucleons may deform or excite as they interact, and they can be virtually 
whatever the dynamics requires, for example, A(or A). This means that we need to consider 
the state mixing effect. In the case of A^A^ interaction, we saw that states mixing plays an 
important role in obtaining the phenomenologically correct potential. We expect the state 
mixing effect to be very important in the A^A^ interaction as well. The state mixing comes 
into effect at the distance where the product ansatz makes no longer sense, so our results 
at short and intermediate distances should be suggestive, although we include state mixing. 
As a guide, we study the effects of the intermediate states A^A, AA^ and AA perturbatively, 
then to second order, the A^A^ interaction is given by 

V(R) = (NN\V{^,C)\NN)+Y. ' <^^l'^(I^-g)l'><-l^(I^-g)l^^) (35) 

Here E^j^ is the two nucleon energy and Eg is the energy of the relevant excited state. 
The first term on the right is the direct nucleon- antinucleon projection of ^(R, C) and it is 
exactly the expression V*^'^''(R). The second term is the correction due to rotational or ex- 
cited states. It is clear from the energy denominator that the second term is attractive. We 
need to evaluate the three sets of matrix elements (A^iV|\/(R, C)|A^A)(A^A|\/(R, C)|A^iV), 
{NN\V{K,C)\AN){AN\V{K,C)\NN) and (A^iV|V(R, C)|AA)(AAJJ/(R,^C)|A^iV), and 



the final result for the first order correction to the A^A^ interaction is 
l^i¥(R) = -^{[\qIP^, + (^P^ + l,Ql)Pl][vl{R) + 2vlm 



19| 



5 ^^3""^ ^27 " 27 

+ {<r' ■ ^')[-Y^QlPo + {f^PN - ^Ql)Pl][vl{R) - vim + (3^^ ■ R^^ ■ R - ■ CT^) 

1 1 fi ^ 

x[-Y^Q^Po + (^P^ - -^Ql)Pl][vl{R)-v,{R)v,m} (36) 
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Here Qn is another finite Nq correction factor Qat = ^(1 — l/Nc){l + h/Nc)- 5 is tlie 
N — A energy difference, which is about 300 MeV, and Pf{T = 0, 1) is a projection operator 
onto the isospin T, PJ" = — Ti ■ T2), P[ = |(3 + n ■ T2). 



B. RESULTS 

For each total isospin T = 0, 1 we parameterize the NN interaction potential by: 

Vnn = yl + ■ + ^^i^^^' ■ R^^ ■ R - ■ (37) 

We now calculate V]^ , Vj , for each isospin T(T=0,1) following the methods outlined 
above. Such a calculation requires considerable computing time. We would like to compare 
the Skyrmion model potentials with the realistic nucleon-antinucleon interaction potentials. 
However, we can not relate our results to the modern interaction potential such as the 



Paris potential 



331] and the Julich potential 3J], since their central parts contain explicit 



momentum dependent terms. For that reason we compare our results with the phenomeno- 
logical potentials of Bryan-Phillips 35| and of the Nijmegen group {36]. These potentials 



provide successful descriptions of both the A^A^ scattering experiments data and the spec- 
trum of resonances, and they are not qualitatively different from each other. At large dis- 
tance all these potentials can be correctly described by the one-boson exchange mechanism 
and the A^A^ potential can be obtained by G-parity transformation of the corresponding 
parts of the A^A^ interaction potential. Using equation of motion and the asymptotic form 
Eq. ([6]) of the chiral angle F{r\ the A^A^ interaction based on the Mod-SKM tends to one 
pion exchange potential in the large distance region fl?! . 

1 „ p—m-^r 

^ ' 47r^2Mjv r ^ ^ ^ ' 

The parameters e,FT^,m are properly chose to guarantee that the long distance tail of the 
A^A^ interaction will agree with the phenomenology. In order to model the annihilation effect 
at short distance, various cut off has been used in the Bryan-Phillips, Nijmegen, and other 
similar potentials. At short distance, the interaction is dominated by the strong absorptive 
potential of order 1 GeV, and it is significantly different from the meson exchange potential. 
Furthermore, the Skyrme model at short distance is no longer meaningful, so we should not 
take seriously the comparison of our results with the phenomenological potentials at 1 fm 
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and less, however the results are still indicative at short distance. We find that the principal 

feature of the phenomenological A^A^ interaction emerges from the careful calculation of that 

interaction based on the Mod-SKM, i.e., the strong central attraction. 

Fig. 2 and Fig. 3 show the central potential calculated from Eq. ( l35i) and the first 

term of the right hand of Eq. (l35l) only. In order to keep the figures clear, we plot the 

potential curves of T = and T = 1 separately. For the case with the nucleon only, the 

results of are independent of the isospin T, and less attractive than the phenomenological 

potentials. When the perturbation corrections due to the effects of the intermediate states 

NA, AN and AA (i.e., A mixing effects) are taken into account, the results of show 

significant attraction effects explicitly, and are closer to the Bryan-Phillips and Nijmegen 

potential. These perturbation results are rather realistic. The effects of A(A) mixing are 

so striking in the case of T = 1 that the perturbation result is more attractive than the 

phenomenological potential for T = 1. Furthermore, we would like to mention that, due to 

isospin conservation, the A^A transition is missing in the T = channel which differentiates 

then the effect of the perturbation result between the T = and the T = 1 channels. As a 

cross-check of our numerical calculation, we reproduce the results of Ref. 19l, the nucleon 
□ 

only results of Ref. [19|] are also shown in Fig. 2 and Fig. 3 in order to illustrate the role of 
B'^B^ term. From these figures it can be seen that the central potentials from the Mod-SKM 
are in better agreement with the phenomenology potentials. 

In Fig. 4 and Fig. 5, we show the T = and T = 1 spin-dependent potentials. In 
these cases the nucleon only potential and the perturbative results are quite similar. From 
1 fm to about 1.5 fm, the potentials from the modified Skyrme model are not so close 
to the phenomenological potentials. Especially in the T = case, both the nucleon only 
and perturbative analysis give a positive spin-spin potential, in contrast to the negative 
values of the phenomenological potentials. It is important to see if the more complete 
Skyrme calculations can repair this disagreement. However, the smallness of the potential is 
reproduced. In our calculation, that small value arises from the cancelations of large terms. 

Fig. 6 and Fig. 7 show the tensor potential V^'". Being similar to case of the spin- 
dependent potential, the nucleon only potential and the perturbative results are also quite 
similar. Particularly at large distance, these results agree with the phenomenological po- 
tential, but the agreement is not so good at intermediate distance. However, the difference 
between the theoretical and the phenomenological results is of the order of 10 MeV, com- 
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pared to the static soliton mass or the nucleon mass which is about IGeV, the difference is 
small enough. Here again, an improved Skyrme model dynamical calculation, going beyond 
the product ansatz, using diagonalization for state mixing and including explicitly the vec- 
tor meson (p, u) and some high derivative terms in the Lagrangian, might lead to a better 
agreement. 



IV. CONCLUSION AND DISCUSSION 

We have shown that the modified Skyrme model with product ansatz can give NN in- 
teraction which is in qualitatively agreement with the phenomenological potential, and it 
provides a better description of the static properties of single baryon than the minimal ver- 
sion Skyrme model. We see that the configuration mixing is very important to be included, 
and we roughly estimate this effect by perturbation theory. More sophisticated method of 
considering the state mixing effect is the Born-Oppenheimer approximation. The poten- 
tial curves in the Born-Oppenheimer approximation are similar to the perturbative results. 



especially for the spin-dependent and the tensor potential |l6l. Il9|. 

To go from this work to a theory that can be confronted with experiment in detail is 
a difficult challenge, i.e, predicting the nucleon-antinucleon scattering cross section, the 
polarization, and the spectrums of the nucleon-antinucleon system etc. There are non- 
adiabatic effects that are particularly important at small R, and there are other mesons 
which should be included in the Skyrme lagrangian. The effects due to vector mesons 
may be particularly important at small distances. Obtaining the static nucleon-antinucleon 
interaction from Skyrme model based on large Nq QCD can be a promising approach. 
We expect we can further discuss whether or not there exists nucleon-antinucleon bound 
state(baryonium) in this framework. 
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FIG. 1: Various terms of the Skyrmion-antiSkyrmion potential potential Eq. (j29p . vi{R) is shown 
by the solid line, V2{R) by the dashed line, V3{R) by the dotted line, V4{R) by the dash-dotted 
line, f5(i?) by the dash-dot-dotted line, and vq{R) by the short dash-dotted line. 
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FIG. 2: The central potential as a function of distance R for the T = channel. The solid 
line is the nucleon only result, the dashed line is the result of the states mixing by perturbation 
theory. The dash-dot-dotted line is the nucleon only potential in conventional Skyrme model [3]. 
The phenomenological potentials based on meson exchange are shown by dash-dotted line for the 
Bryan-Phillips potential and by dotted line for the Nijmegen potential. 
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FIG. 3: The central potential V^F as a function of distance R for the T = 1 channel, labeling of 
the curves is the same as that in Fig. 2. 
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FIG. 4: The spin-dependent potential as a function of distance R for the T = channel, 
labeling of the curves is the same as that in Fig. 2. 
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FIG. 5: The spin-dependent potential, same as Fig. 4 but for T = 1. 
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FIG. 6: The tensor potential Vj^ as a function of distance R for the T = channel, labeling of the 
curves is the same as that in Fig. 2. 
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FIG. 7: The tensor potential, same as Fig. 6 but for T = 1 
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